Damage spreading during domain growth 
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We study damage spreading in models of two-dimensional 
systems undergoing first order phase transitions. We consider 
several models from the same non-conserved order parameter 
universality class, and find unexpected differences between 
them. An exact solution of the Ohta-Jasnow-Kawasaki model 
yields the damage growth law D ~ t*, where <j> — t d ^ 4 in d 
dimensions. In contrast, time-dependent Ginzburg-Landau 
simulations and Ising simulations in d = 2 using heat-bath 
dynamics show power-law growth, but with an exponent of 
approximately 0.36, independent of the system sizes studied. 
In marked contrast, Metropolis dynamics shows damage grow- 
ing via cf> ~ 1, although the damage difference grows as t 0A . 
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The time dependence of fluctuations and correlations 
is a central aspect of the study of equilibrium and 
non-equilibrium phenomena. Fluctuations in the time- 
evolution of a dynamical system have two different ori- 
gins, arising either from the random assignment of the 
initial conditions or from random influences during the 
dynamics - i.e. the (thermal) noise. A particularly in- 
teresting method for isolating the influence of initial con- 
ditions on the subsequent evolution is damage spreading 
0. Consider a kinetic Ising model with site spins ±1, 
evolving via Monte Carlo dynamics. We introduce dam- 
age as follows. Starting with two lattices with identical 
spin configurations we take one system, select at random 
a single site and flip the spin. This introduces a micro- 
scopic fluctuation into the 'damaged' system relative to 
the undamaged one. Both lattices, including the dam- 
aged site, are then evolved according to the same Monte 
Carlo dynamics, by using the same sequence of random 
numbers for both. Any subsequent difference between 
the two is due to the difference in the initial conditions, 
since without the initial damage the two systems are al- 
ways identical. 

One common measure of damage between two Ising 
systems is the Hamming distance between the configura- 
tions, given by 



L d is the number of spins in a hypercubic lattice of size 
L. Another important measure is the damage difference, 
defined by 



D(L,t) = l^\Si(t)-SP(t)\ 



(1) 



r(L,t) 



lj2(So(0)-S D (0)) (Si(t)-SP(t)) 



(2) 



where Si is the state of spin i in the undamaged lattice, 
SP is the state of site i in the damaged lattice and = 



This measures the difference between the amount of dam- 
age that is the same as that initially imposed at site o 
and time t = and the amount that is of opposite orien- 
tation. The term (S o (0) - SP(0)) defines the 'parity' of 
the initial damage at site o and time t = 0. We say that 
damage is positive if it is the same as that initially im- 
posed and negative otherwise. At equilibrium T(t) can be 
related to thermodynamic quantities such as the suscep- 
tibility and correlation functions, provided the damaged 
site is frozen for all time in its t — configuration 
There are no known relationships between damage and 
correlations when the damaged site is allowed to evolve, 
which is the method considered here. Lastly we remark 
that damage evolution in systems far from equilibrium 
has been little studied, and that the relationship between 
damage growth and traditional dynamical properties in 
these cases is unknown. 

In this paper we examine this latter situation, namely 
the evolution of damage in non-equilibrium systems, 
specifically during domain growth. Following a quench 
from a disordered state to a low-temperature state where 
phases coexist, a system dynamically evolves by the for- 
mation and subsequent growth of domains of ordered 
phase. The average size of those domains grows as 
R ~ t n , for late times where n is the growth expo- 
nent, and the morphological evolution involves scale in- 
variance, where all time dependence enters through the 
characteristic diverging length R. An important issue in 
this area is universality: what are the common features 
which characterize a universality class, like n, and how 
are universality classes determined? The most well es- 
tablished universality class in domain growth is model 
A Q , where the scalar order parameter is nonconserved, 
and |j] n = 1/2. Systems in this class include binary al- 
loys undergoing order-disorder transitions, as well as the 
ferromagnetic spin-flip Ising model, and continuum mod- 
els such as the time-dependent Ginzburg-Landau equa- 
tion (TDGL) M and the Ohta, Jasnow and Kawasaki 
(OJK) model km. Herein, we investigate damage evolu- 
tion during domain growth in these latter three models. 



1 



We find some similarities, but important differences, for 
the different models. Hence we suggest that universality 
in domain growth does not include damage evolution. 

We first generalize the concept of damage to allow for 
non-discrete variables. With this definition we are able 
to show that OJK dynamics leads to D ~ t^, with a dam- 
age exponent <f> = d/A in d dimensions. Two-dimensional 
simulation studies of the TDGL equations and of the 
Ising model using heat-bath dynamics both yield power- 
law growth for damage. However, <\> w 0.36, significantly 
smaller than the OJK model. More remarkably, Ising 
simulations using Metropolis dynamics in d = 2 show 
<f> ~ 1. It should be noted that in all the models studied 
damage disappears with high probability at long times. 
But for some realizations of the initial conditions and 
thermal noises, it attains macroscopic size. The power 
laws above apply to the damage averaged over samples 
of the full ensemble of realizations. 

The OJK model || was introduced as an approxima- 
tion to the Allen-Cahn equation describing interface mo- 
tion during phase ordering with a non-conserved order 
parameter. It is considered an important model in its 
own right H as it readily yields analytic predictions for 
quantities such as the growth exponent and scaling func- 
tions. 

Consider the evolution of a system described by an 
order parameter ^(x, t) which takes the value ip eq in one 
of the phases and — ip eq in the other. In the OJK model 
the evolution of ^(x, t) is given by 

u(x t) 

V>(x, t) = -0 eg sign(u(x, t)) = ip eq ' ( , (3) 

where the auxiliary field w(x, t) satisfies the diffusion 
equation: u(x, t) — V 2 u(x, t). Without loss of generality 
the values of ip eq and of the diffusion coefficient are taken 
here to be unity. The initial condition Uo(x) = u(x, t — 
0) is usually taken as a random Gaussian field of zero av- 
erage and short-range correlations. The exact solution of 
the diffusion equation is it(x, t) = J <fxG(x— x', t)uo(x'), 
where G(x ,t) is a Green function. Substitution of this 
expression into Eq. (^) leads to an equation of evolution 
for tp(x, i) representing the growth of domains of typical 
size R(t) ~ t 1 / 2 §. 

We need continuum definitions for the damage measure 
D(L,t) and for the initial damage consistent with those 
used for the Ising model. We define the initial damage 
in terms of the initial unperturbed field Uo(x). Several 
definitions of the 'damaged' field Uq are possible: we 
found all to give essentially identical results. We present 
detailed results for the case in which Uq(x) = — uo(x) 
in a damaged region of volume Vd and Uq'(x) = ito(x) 
outside. We define the damage measure by 

D{L,t) = \ f dx |V>(x,t)-V^(x,i)| (4) 
= U Ld ~ j[/x^(x,i)/(x,t)). (5) 



The calculation of the average damage (ensemble- 
averaged over different realizations of the random initial 
conditions) then reduces to the calculation of a correla- 
tion function and a subsequent integration. The detailed 
calculation will be presented elsewhere: here we quote 
the main results. 

In the regime of long times and large L and with the 
initially damaged volume much smaller than system size 
we find the following scaling result 

(D(L,t)) — F(t/L 2 ), (6) 
where the scaling function F is 

F (r) = 2V J° J— f dx G(x,r). (7) 

I d is the the unit hypercube. The crucial and disputable 
element here is the L d l 2 factor. This arises from the 
assumption of Gaussian de-coupling in the OJK model. 
The scaling form in the limit t <C L 2 is 

(D{t^L 2 )) ~ti. (8) 

Plots of (D(L,t)) from Eqs. (0) and @j for d = 2 are 
shown in scaled form in Fig. [if The power-law regime 
is clearly seen, as are the corrections to scaling at short 
times. 

We also numerically iterated the OJK equations on a 
lattice to visualize the evolution of the damage. We found 
that damage, when it exists, evolves along the interface 
between the two phases. Damage then either disappears 
as the interfaces disappear, or it grows to encompass the 
entire system. In addition we confirmed that only one 
parity of damage dominates: in particular if we damage 
only a single lattice site the propagated damage is always 
of the same parity as initially imposed. This result can be 
verified analytically. Therefore in OJK, with this dam- 
age prescription, the damage difference and Hamming 
distance are identical. 

The OJK damage definition (Eq. (|J)) also applies to 
the TDGL case. The initial damaged configuration in 
this case was chosen to be ip D (x, t = 0) = —ip(x,t = 0) 
inside a small damaged region and^x,*) = ip(x,t = 0) 
outside, where again ^(x, t — 0) is a Gaussian random 
field with short range correlations. Both ip and ip D are 
evolved according to the TDGL equation: 

V>(x, t) = V>(x, t) + VV(x, t) - <A(x, t) 3 . (9) 

Addition of a noise term does not significatively change 
the results. The average damage evolution is shown in 
Fig. ^|. The power-law damage growth regime is clearly 
present and persists up until the system magnetization 
begins to saturate. We obtain <f> = 0.36 ± 0.02, which 
is smaller than the OJK prediction of d/4 = 1/2. The 
errors are estimated by one standard deviation of the 
statistics of the data. 
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As with the OJK model, we found that the TDGL 
equations, solved on a discrete lattice, only propagate 
damage of the same parity as that initially imposed, pro- 
vided we only damage a single lattice site. If we damage 
a larger domain there is transient regime during which 
both parities of damage compete, but eventually one dis- 
appears and damage follows <j> » 0.36. 

Ising simulations, using both heat-bath and Metropo- 
lis dynamics, were performed on two-dimensional square 
lattices, using a multi-spin coding algorithm described 
previously |||| . However the block-spin coding trick ||] 
was avoided, as some kinds of parallel updating have been 
found to affect damage dynamics |I^,|llJ. Spins to be 
evolved were chosen at random from amongst all spins, 
a single time step (mcs/s) being L 2 such attempts. The 
initial state was an infinite temperature equilibrium con- 
figuration. 

Typical results for damage evolution via heat-bath dy- 
namics are shown in Fig. |2|, for two system sizes and 
for temperatures T = 0.8T C and T = 0. The early-time 
transient is very short, and is quickly followed by a slow 
power-law growth regime. The power-law exponents are 
0.36 ±0.02, identical with the TDGL results and smaller 
than the OJK prediction of 1/2. It also appears that the 
system size does not strongly affect this exponent: sys- 
tems of size L = 48 have essentially the same slope as 
those four times as big. There is no evidence for faster 
growth at late times: the growth simply slows as finite- 
size effects are encountered, and then saturates. Temper- 
ature also does not appear to affect the exponent in the 
power-law regime. In general increasing the temperature 
decreases the total damage. This is consistent with pre- 
vious studies of damage in equilibrium systems using the 
heat-bath algorithm. 

We cannot rule out the possibility that the 0.36 expo- 
nent is only effective, and that systematic errors preclude 
our observation of d/A. Our estimated errors are statisti- 
cal, and comparable systematic errors are possible. How- 
ever, we have made a careful study of systems of different 
sizes, and our numerical results for (j) are not consistent 
with the exact result of the OJK model of d/A. 

It was previously demonstrated that the heat-bath al- 
gorithm can only propagate damage of the same parity 
as that initially imposed 0. Therefore as in the pre- 
vious two cases the Hamming measure and the damage 
difference are identical. 

We obtain remarkably different behavior using the 
Metropolis algorithm. Typical Metropolis dynamics re- 
sults are shown in Fig. ^, where again we have plotted re- 
sults demonstrating the effect of varying system size and 
temperature. The results are different from all results 
presented above. Even at earliest times the growth of 
damage is faster than i 1 / 2 , with this exponent increasing 
with time. At intermediate times this damage saturates 
at <f> ~ 1, i.e. approximately twice as fast as seen by the 
other dynamics. Simulations to much later times indi- 
cate no deviation from this dependence until finite-size 
effects are encountered and growth stops. 



Equivalence between the three numerical models can 
be recovered if we consider the damage difference, Eq. 
(||). A typical plot of the damage difference is shown 
in Fig. H for systems of size L — 48. Although quite 
noisy, these data clearly show a power-law growth slower 
than t 1 / 2 . A best fit yields 0.40 ± 0.04, consistent with 
the TDGL and Ising heat-bath results. Of course we 
again cannot rule out that the possibility that this expo- 
nent approaches 0.5 for larger systems, indicating strong 
finite-size effects. However, we see no evidence for such 
systematic errors in the range of sizes studied (L = 48 to 
L = 192). 

The reason for the surprising differences between the 
Ising simulations lies, at the microscopic level, in the dif- 
ference between the two dynamics. Metropolis dynamics 
allows both parities of damage to exist. By visualizing 
the damage evolution we find that Metropolis dynamics 
produces significant amounts of both damage. Individ- 
ually averaging these two parities of damage shows that 
both are dominated by <f> ~ 1. Damage of the same par- 
ity as the initial damage shows a large initial transient 
related to the early-time growth of the damage differ- 
ence, before crossing over to the faster growth, while the 
opposite-parity damage quickly assumes <f> ~ 1. 

In conclusion, we have presented numerical evidence 
that non-equilibrium damage spreading lacks the univer- 
sality of domain growth. In particular we have shown 
that the TDGL model and the Ising models are equiva- 
lent, but only when the damage difference is used as the 
measure. More importantly our work suggests that the 
damage algorithm breaks the universality between the 
TDGL, Ising, and OJK models. More detailed compar- 
isons between these models, including the consideration 
of finite-size effects, finite-size scaling, and the full deriva- 
tion of the OJK analytical results will be presented in a 
future publication. 
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FIG. 1. The d — 2 OJK prediction for a initially damaged 
region of Vd = 4 and three different system sizes L. The 
scaling region and the damage growth via are evident, as 
are corrections to scaling at short times. 

FIG. 2. Damage evolution for a 2d quenched Ising model 
using heat-bath dynamics, showing the effect of temperature 
and of system size. The growth of damage is evidently alge- 
braic, even for the smaller systems (L = 48). The exponents 
are 0.36 ± 0.01 for L = 48, T = 0.8T C (24576 runs averaged); 
0.37 ± 0.03 for L = 192, T = 0.0T C (3072 runs averaged) and 
0.36 ±0.01 for L = 192, T = 0.8T C (3072 runs averaged). The 
dotted line shows the TDGL damage evolution averaged over 
6400 realizations of the initial conditions, arbitrarily scaled 
to fit on the same scale as the Ising data. The system size is 
L — 90.5 and the initially damaged volume is Vd = 0.5. Note 
the discrepancy with the t d ^ A = t 1 ^ 2 OJK prediction. 

FIG. 3. Damage evolution for a 2d quenched Ising model 
using Metropolis dynamics, showing the effect of temperature 
and of system size. The regimes for damage are only evident 
for the larger system, and at late times. The solid line shows 
a damage evolution oc t. The damage difference is also plot- 
ted, and the slow-growth behavior is evident. The late-time 
tangent to the growth gives an exponent of <f> — 0.40 ± 0.04. 
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